Graphical Method for Effective Interaction with a New Vertex Function 
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Introducing a new vertex function, Z{E), of an energy variable E, we derive a new equation for 
tlie effective interaction. Tiie equation is obtained by replacing the Q-box in the Krenciglowa-Kuo 
(KK) method by Z{E). This new approach can be viewed as an extension of the KK method. We 
show that this equation can be solved both in iterative and non-iterative ways. We observe that the 
iteration procedure with Z{E) brings about fast convergence compared to the usual KK method. It 
is shown that, as in the KK approach, the procedure of calculating the effective interaction can be 
reduced to determining the true eigenvalues of the original Hamiltonian H and they can be obtained 
as the positions of intersections of graphs generated from Z{E). We find that this graphical method 
yields always precise results and reproduces any of the true eigenvalues of H. The calculation in the 
present approach can be made regardless of overlaps with the model space and energy differences 
between unperturbed energies and the eigenvalues of H. We find also that Z{E) is a well-behaved 
function of E and has no singularity. These characteristics of the present approach ensure stability 
in actual calculations and would be helpful to resolve some difficulties due to the presence of poles 
in the Q-box. Performing test calculations, we verify numerically theoretical predictions made in 
the present approach. 

PACS numbers: 21.60.De, 24.10.Cn, 02.30.Vv, 02.60.Cb 



I. INTRODUCTION 



In nuclear, atomic, and chemical physics, it is often 
useful to introduce an effective interaction acting in a 
chosen model space (P space) of tractable dimension. In 
nuclear physics, much effort has been made both as re- 
gards formal theories and their applications [H-Q- R-e- 
cently the effective interaction method has been applied 
to new fields of many-body physics, such as quantum 
dots [§, and many-boson systems [lH . 

Among many approaches, we here direct our atten- 
tion to the Krenciglowa-Kuo (KK) [l^, and the Lee- 
Suzuki (LS) 14, 15] methods. These two methods are 
constructed in terms of the so called Q-box as a building 
block of formulation. The KK approach has very sim- 
ple structure and the effective interaction is obtained in 
an iterative way. If the iteration converges, in almost 
all numerical calculations, eigenvalues are given for the 
states which have the largest overlaps with the chosen 
model space. On the other hand, the LS method repro- 
duces eigenvalues for the states which lie closest to the 
chosen unperturbed energy. Originally the LS method 
had been presented to resolve the difficulty of divergence 
in the perturbation expansion. The LS method is rather 
complicated in structure and higher derivatives of the Q- 
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box with respect to starting energy are necessary if one 
wishes to obtain more accurate solutions. 

Both of the two theories yield only certain of the true 
eigenvalues of the original Hamiltonian. This restriction 
is not desirable. In a formal point of view, the Q-hox it- 
self contains information regarding all of the true eigen- 
values. For a given model space of dimension d, there 
would be a method of reproducing any d eigenvalues 
among all the true eigenvalues. 

In many cases of applying the effective interaction the- 
ories an iteration or a recursion method has often been 
employed. In general, the convergence in such a method 
depends strongly on the properties of the eigenstates of 
the Hamiltonian H. In actual calculations the infor- 
mation on the true eigenstates is not given beforehand. 
Therefore, it is impossible to control the convergence in 
the iteration. In many cases we cannot know whether the 
iteration is convergent or not before starting calculations. 
Even if the iteration converges, we do not know which 
eigenvalues are reproduced in all of the eigenvalues. An- 
other difficulty encountered in actual calculations is the 
pole problem. The Q-box itself has poles at the energies 
which are the eigenvalues of H in the complement space 
(Q space) . The presence of poles causes often instability 
in numerical calculations. 

We shall show that it is indeed possible to resolve these 
difficulties by introducing a new vertex function, Z{E), 
in place of the Q-box by Kuo et al. [H, [l^l • Preliminary 
version of the present work has been reported in Ref . [iq . 
Very recently Dong, Kuo, and Holt have followed the 
present approach and applied to the actual calculations 
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of the shell-model effective interactions flT'l . These works 
have shown that the present method has a possibility of 
providing a suitable framework with Z{E) as well as the 
KK approach. 

The organization of the present paper is as follows: 
In Sec. [n] we outline the Q-box formalism and the KK 
method. A new vertex function operator, Z{E), is in- 
troduced and some of its mathematical properties are 
clarified. The algorithm of the calculation procedure is 
given by applying the secant and the Newton-Raphson 
methods. In Sec. IIIII test calculations are made in or- 
der to assess the present approach. We examine whether 
the theoretical predictions are verified or not in an ex- 
actly solvable model. Concluding remarks are given in 
SecHVl 



II. FORMULATION 

We consider a general quantum system which is de- 
scribed by a Hamiltonian H. We write an eigenvalue 
equation with the eigenvalues {Ek} and the eigenstates 
{|*fe)} as 



fc = 1,2,. 



(1) 



The Hamiltonian H is supposed to be composed of the 
unperturbed Hamiltonian Hq and the perturbation V, 
i.e., H = Hq + V. We decompose the entire Hilbert space 
into the model space (P space) and its complement (Q 
space) with the projection operators P and Q, respec- 
tively. Basic properties of the projection operators are 
P + Q ^ 1 and PQ = QP = 0. 

We here assume that Hq is decoupled between the P 
and Q spaces as 



Ho = PHoP + QHoQ. 



(2) 



In the present work we consider a case that the P-space 
states have a degenerate unperturbed energy Eq, i.e.. 



PHoP = EoP 



(3) 



Then the P-space eigenvalue equation is written with the 
effective interaction R and the eigenstate \(j)k) as 



{EoP + R)\cf>k) ^ EM 



fc = l,2,...,d, (4) 



where d is the dimension of the P space. In the above 
d eigenvalues Ek 's should agree with d of the true eigen- 
values in Eq. ([T]). Various solutions for R are possible, 
and many theoretical frameworks have been given for ob- 
taining R. Probably the most widely applied effective 
interaction is given by imposing the condition that the 
model-space eigenstate \(j)k) in Eq. Q should agree with 
the P-space component of the true eigenstate \^>k) of H, 
i.e., \(j)k) = P\'^k)- This restriction on \(f>k) leads to the 
standard non-Hermitian form of R [ISj] . 



Solutions for effective interaction in the Q-box 
formalism 



Among many approaches to the effective interaction R 
we here discuss the KK formalism [l3,[l3|- Originally the 
KK method is based on a diagrammatic representation of 
the effective interaction, which has been known as the Q- 
box folded-diagram method originated by Kuo et al. [l9| . 

For obtaining R, one first calculates the vertex function 
called the Q-box which is defined as the sum of all the 
linked and non-folded diagrams. Next one should add the 
folded diagrams, which can be carried out rather simply 
by applying the energy-derivative expression of the Q- 
box J_9,]- Originally the KK approach was proposed to 
derive an effective interaction acting among a few valence 
particles outside the closed-shell core. However, in the 
present work, we consider a general quantum system and 
wish to reproduce d total energies of H by introducing 
an effective interaction. For this case the KK method 
can also be applied by defining the Q-hox in an operator 
form as 



Q{E) = PVP + PVQ 



1 



E-QHQ 



QVP, 



(5) 



which is a function of an energy variable E. Equation ([5]) 
is equivalent to the energy-dependent form of the effective 
interaction given by Bloch and Horowitz in the many- 
body perturbation theory [20| . 

With the Q-box the effective interaction R can be ex- 
panded into 

R = Q + QiQ + QiQiQ + Q2QQ + ---, (6) 
where Q = Q{Eq) and Qm = Qm{Eo) with 

Here Eq is the starting energy or the degenerate unper- 
turbed energy as given in Eq. (jS]). 

If the series expansion in Eq. ^ converges, R is given 
in a formal way by 



R = J2QiEkMk){^k\, 



(8) 



k=l 



where E^ and |0fc) are given in Eq. and {(j>k\ is 
the biorthogonal state defined through the orthogonality 
{4'k\4>k'} = 5kk' ■ We here note that, as seen in Eqs. (|4]) 
and ([H]), the derivation of the effective interaction R is 
equivalent to determining d eigenvalues {Ek\ and the 
corresponding P-space eigenstates { I (^fc)}. Using Eqs. (g]) 
and dl]), the effective interaction R is rewritten simply as 



R^Y.^Ek~Eo)\cpk){k\ 



(9) 



fc=i 
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and the model-space eigenvalue equation with the effec- 
tive interaction R is expressed as 



[EoP + Q{Ek)]\<l>k) ^ EM 



(10) 



From Eqs. ([8]) and ^TU\\ we understand that R can be 
given by calculating the Q-box at the true eigenvalues 
{Ek} of H as starting energies which are determined self- 
consistently. 

We here show that there is another way of deriving the 
effective interaction R. We first note that the Q-box is a 
function of E and resultantly the eigenvalues of EqP + 
Q{E) are also functions of E. We write the eigenvalue 
equations in the P space for an arbitrary energy E as 

[E^P + Q{E)]\^„,) - G,n{E)\ilj„,),m = 1, 2, • • • , d. (11) 

Since the P space is d dimensional, we have d eigenvalues 
denoted by {Gm{E)] ni — 1, 2, • • • , d}, which we label in 
order of energy as Gi{E) < Ga (£;)<••• < Gd{E). It 
may be clear from Eqs. ([TU]) and (ITT|) that the eigenvalues 
{Ek} in Eq. ([T0| can be given by solving the following 
equations 



G^iE)^E, 



1,2, 



(12) 



Various mathematical methods have been known to solve 
such equations. 

We should note that a set of the equations are 
independent of the properties of the eigenstates of H, 
such as P-space overlaps or energy spacings. Therefore, 
in principle, it is possible to reproduce the true eigenval- 
ues {Ek} of H more than in the usual KK method based 
on Eq. pUj) . However, as seen in Eq. ([5]), there appear 
poles in the Q-hox when E approaches one of the eigen- 
values of QHQ. The poles in Q{E) induce also the poles 
in Gm{E) in Eq. (Ill|) . Such a situation causes instability 
in numerically solving Eq. for {Ek} around the pole 
positions. 



B. Extension of the Krenciglowa-Kuo method 

In order to resolve the pole problem we introduce a 
new vertex function of an energy variable E, which is a 
P-space operator defined in terms of the Q-box and its 
energy derivative as [ll| 



Z{E) 



1 



l-Qi{E) 



Q{E) - Qi{E){E ^ Eo)P . (13) 



Hereafter we shall refer to Z{E) as the Z-hox. We note 
here that the Z-hox agrees, at E = Eq, with the first- 
order recursive solution in the LS method UM ■ In a recent 



paper of Dong, Kuo, and Holt [17[ the definition of Z{E) 
is given in a more general case with a non-degenerate P- 
space unperturbed Hamiltonian PHqP. In this case one 
should replace EqP by PHqP. As in Eq. ((TTl) with the 
Q-box, we consider an eigenvalue problem 

[EaP + Z{E)]\^m) ^ Frr,iE)\^^^), m = 1, 2, • • • , d,(14) 



where {Fm{E); m = 1,2, ■ ■ ■ , d} are d eigenvalues which 
are functions of E. We here label {Fm(£');m = 
1,2, - ■■ ,d} in order of energy as Fi{E) < F2{E) < ■ ■ ■ < 
Fd{E). ^ 

The Z-box and the associated functions {F,n{E)} have 
the following properties: 

(i) Using Eqs. ©, dH), and ([13]), we have, for the 
P-space eigenstates {|0fc)} in Eq. (fTO|) . 

d 

J2z{EkMk){k\ 



k=l 



(Ek 



R-QiiEk)-R \^k){M 



R. 



(15) 



The above fact means that, replacing Q{E) by 
Z{E) in Eq. a new solution for the effective 
interaction R can be derived as 



REKK = Y.^(^kMk){4>k 



(16) 



k=l 



and equivalently Eq. pU)) with the Q-box is re- 
placed by 



[EoP + Z{Ek)]\(f>k) 



Ek 



k\Vk) 



(17) 



From the above relations between two approaches 
with Q{E) and Z{E) we may call Rekk hr Eq. 
the extended Krenciglowa-Kuo (EKK) solution. In 
the same way as in Eq. p^ . the true eigenvalues 
{Ek} can be given by solving the equations 



F^(E)^E, TO = 1,2, 



(18) 



(u) Using Eqs. (|5l). (fT)). and p3| we can derive a formal 
expression for the energy derivative of Z{E) as 



dZ{E) _ 2 
dE ~ l-Qi{E) 

x[Z{E)-{E 



2{E) 

Eo)P]. 



(19) 



In the paper of Dong et al. [17[ , the above expres- 
sion has also been given for a general case with the 
non-degenerate PH^P. If E is one of the true eigen- 
values {Ek} satisfying Eq. pT|) . we see that the en- 
ergy derivative of Z{E) becomes zero at E — Ek, 
namely 



dZ{E) 



dE 



\<f>k)^0, k = l,2,---,d. (20) 



E=E., 



Resultantly we have for the energy derivative of 

F,niE) 



dF,r,{E) 



dE 



= 0, k = l,2,---,d. 



(21) 



E=Ek 
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These results for the energy derivatives have been 
pointed out in the previous paper [T^. 

(iii) We discuss here some problems associated with the 
poles of Q{E). First we consider the eigenvalue 
equation for the Q-space Hamiltonian QHQ writ- 
ten as 



QHQ\q)=e,\q), 



(22) 



where Sg and \q) are the eigenvalue and the eigen- 
state, respectively. It may be clear from Eq. ([5]) 
that Q{E) has a pole at E — Sq. We define a P- 
space operator Xq with the Q-space eigenstate \q) 
in Eq. as 



Xq = PV\q){q\VP. 



(23) 



We write the eigenvalue equation for Xq with an 
eigenvalue 



Xqlfi) = Xf,\fl). 



(24) 



From the definition of Xq in Eq. ([23]) we see that 
Xq is Hermitian and positive semi-definite, that is, 
Xq has positive or zero eigenvalue because of the in- 
equality = > 0. We can further prove 
that there is only one eigenstate, denoted by |/io), 
with a positive eigenvalue and that the eigenvalues 
of all the other eigenstates are zero. The only one 
positive eigenvalue of Eq. is given by 



= J2\(p^\^\i)\'' 



(25) 



i=l 



where {\pi); i — 1,2, ■ ■ ■ , d} are the basis-state vec- 
tors of the P space. The proof is as follows: We 
condsider a matrix representation of Xq with the 
matrix element iXq)ij — {pi\V\q){q\V\pj). Then 



we obtain the trace of Xq as 



TrXq = = 



(26) 



1=1 



We note that the state vector |/io) can be written 
explicitly as 

1^^^ ^ [{pi\V\q)AP2\V\q),---APd\V\'l)V ^ (27) 



N 

where the symbol T means transpose of (1 x d) 
matrix, and ^/N is the normalization factor with 
= J2i=i KPil^k)^- This state |/xo) becomes an 
eigenstate of Xq and has a positive eigenvalue as 

= (J2\{p^\V\q)\'^ \f^o) 

= (28) 



The eigenvalue coincides with the trace of Xq. 
Recalling a well-known theorem on the trace of a 
square matrix, we see that the trace of Xq should 
be equal to the sum of the eigenvalues of Xq. It 
follows immediately that all the other eigenvalues 
except Xfj,^-, become zero. 

(iv) Next, we prove that Eq. ([T7|) or (fT8| has additional 
solutions other than the true eigenvalues {Ek}. We 
note, in the vicinity of a pole of Q{E) where E — 
A with a small deviation A, that Q{E) in 



Eq. ([5]) is expressed as 



QiSq + A) 



PVP+^ 



^7 



x„ 



A 



(29) 



Here we have used the definition of Xq in Eq. 
The operations of Q{eq + A),Qi{eq + A), and 
(32 (Eg + A) on the state |/io) yield in leading or- 
der, respectively 



Q2(£q 



A)|/io) - 
A)lAio) - 
A)|/io) - 



' A2 

A3 



Imo)- 



(30) 
(31) 
(32) 



Using Eq. (|T3|) and the above relations, we obtain 

a.t E ^ Eq + A 

[Eo + Z{eq + A)]\fio) 



''MO 



(e, -f A - Eo 



A2 



Imo) 



(33) 



(34) 



A 

= (£, + 2A)|/io). 

Then we have in the limit of -B — !■ e, 

[Eq + Z{Sq)]\lJ.o) = £q\^J■a)■ 

This fact means that the state |/io) and the energy 
Eq can be an additional solution to Eq. ([T7)) . This 
also means that the pole energy Sq satisfies Eq. (|18p . 

(v) We further consider the energy derivatives of Z{E) 
and Fm{E) at the pole positions. Subtraction of 
Eq. ^ from Eq. ^ follows 



dZ{E) 



dE 



|/io> = 2|/io). 



(35) 



E=e„ 



We can derive the above result, in another way, by 
using Eqs. dH]) and ((M)) 



dZ{E) 



dE 



Imo) 



E=e^+A 

2 



l-gi(e, + A) 



Q2{eq + A)■A\^iQ). (36) 
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Then, using Eqs. ([3T|) and (|32)) and taking the hmit 
of A — > 0, we obtain 



dZ{E) 



dE 



= hni 2 



A3 



A|/^o) 



2|Aio). 



As a consequence we have 



dE 



2. 



(37) 



(38) 



This property of the energy derivative of Fm{E) 
at a pole position is quite contrast to that for a 
true eigenvalue in Eq. (PT|) . These properties of 
the energy derivative of Fm{E) given in Eqs. (pij) 
and ([55]) can be used to distinguish the pole energy 
solutions from those of the true eigenvalues of H . 
It is also noted, from Eqs. (IMl) and (pS)) . that the 
function Fjn{E) is a well-behaved function of E and 
has no singularity even at the pole energies. 



C. Iteration methods for the effective interaction 

Iteration methods have been employed quite often for 
obtaining the effective interaction or equivalently for de- 
termining the true eigenvalues oi H. In the previous 
subsections we have given four basic equations, namely 
Eqs. (Uni) and ^ for the KK and Eqs. ^ and ^ for 
the EKK methods. In the KK approach Eqs. pU|) and 
([T2t lead, respectively, to the following iterative equa- 
tions: 



[EoP 



= E, 



(n+l)| ,(n+l) 



and 



(")^ 



k 



1,2, 



\q>r>) (39) 



(40) 



In the same way two iterative equations can be derived 
from Eqs. pTl) and (IT51) for the EKK method, respec- 
tively as 



[EqP + Z{E, 



E 



(n+l). 



and 



= F™(i?(")), m = l,2,...,d. 



) (41) 



(42) 



The iterative equation ([5^ has long been used as a stan- 
dard method in the KK approach. The convergence con- 
dition has been investigated in many of the theoretical 
and numerical studies [l^, [H, HH . From these studies it 
has been known that, if the iteration converges, the KK 
approach in Eq. pop derives d eigenvalues of the eigen- 
states of H with the largest P-space overlaps. However, 




FIG. 1: Graphical representation of the modified secant 
method. The Li is the straight line passing through the points 
A and B. The L2 and L3 are the tangents to the function f{E) 
at the points A and B, respectively. See text for detailed ex- 
planation. 



it has also been known that the iteration in Eg. ([39)) 
does not always converge. Recently Takayanagi [22] has 
pointed out an exceptional case that the KK method 
does not reproduce states which have the largest P-space 
overlaps even if the iteration converges. Although the 
KK iteration method has been applied widely and has 
brought remarkable results in actual calculations, the rig- 
orous condition of convergence for Eq. (15^ has not yet 
been made clear. 

The convergence condition for the iteration in Eq. PO)) 
requires mathematically that the energy derivative 
of Gm{E) at the solution E — Ek should satisfy 
\dGm{E)/dE\E=Ek < 1- In general, this convergence 
condition is considered to be satisfied by certain of the 
true eigenvalues oi H. It implies that this iteration 
method is restrictive for the purpose of reproducing the 
solutions of {Ek} as many as possible. 

The iteration in Eq. (|4T|) leads to a new scheme of the 
calculations of {Ek} and i?EKK in Eq. ([T5)) . As in the case 
of the KK iteration in Eq. the convergence condition 
in Eq. (1411) is quite complicated and has not yet been 
made clear. However, the iteration in Eq. (j4ip would be 
applicable in actual calculations. In the paper of Dong 
et al. 17;], they have employed this iteration scheme to a 
model problem and actual calculations of the shell-model 
effective interaction. They have concluded that the two 
iterations in Eq. ((23) for the KK and Eq. (gT]) for the 
EKK methods are both suitable and efficient. 

Another iterative equations in the EKK method is 
given in Eq. (j42p . The convergence condition for this iter- 
ation is given by \dFm{E) / dE\E=Ek < 1- It is clear from 
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the property of the energy derivative of Fm{E) in Eq. ((2T|) 
that the convergence condition is satisfied for any of the 
true eigenvalues {Ek}. The iteration in Eq. (|42l) con- 
verges surely to one of {Ek}. We also note that this 
iteration in Eq. (|42|) never reaches the pole-energy so- 
lutions due to the property in Eq. (|38l) of the energy 
derivative of Fm{E). We can show that the iteration 
in Eq. (H^ is equivalent to that in the Newton-Raphson 
method which is often used to solve non-linear equations. 
The Newton-Raphson method has been known to derive 
quadratic convergence, that is, the number of correct dig- 
its is doubled at each step of iteration. On the contrary 
the usual methods, including the KK and the LS meth- 
ods, derive single (or linear) convergence. The iteration 
in Eq. (|42p guarantees surely faster convergence than the 
usual iterations given so far. 



Graphical method for effective interaction with 
the Z-box 



D. 



We present a new method for solving Eq. pa)) derived 
on the basis of the Z-box and its associated functions 
{Fm{E)}. The solutions can be obtained by finding the 
energies of the intersections of two graphs, y = Fm{E) 
and y ~ E. As we have already shown, the Z-box has 
no singularities at the pole positions of the Q-box. Thus, 
{F„i{E)} can be considered to be well-behaved functions 
of E. These characteristics of the functions {Fm{E)} 
enable us to employ some of the mathematically well- 
established methods for solving non-linear equations. 

The procedure of solving Eq. (|18|) in the present ap- 
proach, which we call the graphical method, is given 
as follows: We first draw graphs oi y = Fm{E) for 
m = 1, 2, • • • , d and y = E. The energies at the intersec- 
tions of these graphs become the solutions to Fm{E) = E. 
Furthermore, the energy derivatives of {Fm{E)} at the 
intersections distinguish the pole-energy solutions {Sq} 
from those of the true eigenvalues {Ek}- As was proved 
in the previous subsections, the energy derivatives, de- 
noted by {F^{E)), should be F!,^{E) = Q iov E ^ Ek 
and Fl-^{E) = 2 ioi E = Eg. With these simple regula- 
tions for {Fl^{E)} we can easily specify the solutions of 
the true eigenvalues {Ek}. 

We next figure out, from the intersections of the 
graphs, roughly estimated solutions for {Fm{E) = E}. 
Starting with these approximate solutions, we proceed to 
make further precise calculations for the solutions. For 
this purpose we employ a combined method of the secant 
and the Newton-Raphson methods The algorithm 

of this method is as follows: Let us define a function f{E) 
as f{E) = E — F{E), where F{E) is one of the functions 
{Fjn{E)}. The solutions to Fm{E) = E are obtained as 
the roots of f{E) = 0. Suppose that Ea and Ei, be two 
values which bracket one of the roots and satisfy Ea < Ej, 
and f{Ea)f{Ef,) < 0. We suppose also that f{E) be a 
monotone function on the interval Ea < E < Ft,. These 
suppositions on Fa and Ei, follow that only one root ex- 



TABLE I: Convergence of the eigenvalues of the lowest- 
lying and second lowest-lying states obtained in the iterations 
^(n+i) ^ Gm{E^"'>) and E'^"+^^ = for the KK and 

for the EKK methods, respectively, for the model Hamilto- 
nian with the strength x = 0.05. Correct digits in the KK 
and the EKK methods are given for n, namely the number 
of iterations. The starting energies {E[^\ E^^) are taken to 
be (0.0,0.0). The notation c indicates convergence to more 
than fifteen decimal places. The exact eigenvalues here are 
El = 0.8904504858869942 and E2 = 2.2156808150096040. 



Tl 






1 


0.9- ■ ■ 


0.89--- 


2 


0.88- ■ ■ 


0.890450- - - 


3 


0.8904- • • 


0.89045048588699- ■ • 


4 


0.890450- ■ ■ 


c 


5 


0.8904504- ■ ■ 


c 


6 


0.890450485- - - 


c 


7 


0.8904504858- - ■ 


c 


1 


2.2- ■ ■ 


2.2- - - 


2 


2.21- ■• 


2.21568- ■ - 


3 


2.215- - - 


2.2156808150- - - 


4 


2.2156- - - 


c 


5 


2.21568- - ■ 


c 


6 


2.215680- ■ ■ 


c 


7 


2.2156808- ■ ■ 


c 


8 


2.21568081- - - 


c 


9 


2.215680815- - - 


c 



ists on the interval Ea < E < E\j. As shown in Fig. [1] 
we further determine five points according to 



ei 

63 

64 
65 



= Ea^ 

= ElL 

= Ea 

= Eh- 



■ f{E„) ~ Eb - f{Fa) 
f{Et) - f{Ea) ' 
f{Ea) 
f'iEaY 

riEbY 



Eh 



(43) 
(44) 

(45) 

(46) 
(47) 



The value 62 in Eq. (|44l) is the better approximate so- 
lution in the usual secant method. The 63 and 64 are 
the approximate values in the Newton-Raphson method. 
We can easily select two values, and ej, among 
{ei, 62, • • • , 65} such that they are neighboring on the E- 
axis and satisfy Ea < Ci < ej < E^ and f{ei)f{ej) < 0. 
These conditions mean that one of the roots exists on 
the interval Ci < E < ej. We then replace Fa and Ei, 
as Ea = Ci and Ei, = Cj and repeat the procedure again. 
This modified secant method derives surely one of the 
roots on the interval Fa < E < Ft,. We also note that 
this method will lead to fast convergence, because we cal- 
culate approximate solutions according to the Newton- 
Raphson method. 

We here wish to emphasize that this method is not an 
iteration method. The procedure in this method always 
derives a convergent solution to f{F) = on the inter- 



7 



Gm{E) 




FIG. 2: Dependence of {G^(E)} on E with x = 0.2. The 
graphs oi y = Gi (E) and y = G2 (E) are shown in solid and 
broken hnes, respectively. The graph of y = E is also shown. 

val Ea < E < Eh and never reaches any of the solutions 
outside the interval. Therefore, we can always control 
convergence and select solutions to be reproduced by se- 
lecting starting values Ea and Eb properly. Only one 
problem in the present method is how to choose Ea and 
Eb- Then the use of the graphs of {y = Fm{E)} would 
be helpful. If the graphs are drawn accurately, we can 
determine easily these starting energies Ea and Eb- 

III. TEST CALCULATIONS 

In order to assess the present method we consider a 
model problem for which exact results can be obtained 
easily. The model we adopt here [llj is a slightly modi- 
fied version of the one which was introduced many years 
ago to study the intruder state problem [2^, |25|] . The di- 
mensions of the entire space and the P space are four and 
two {d = 2), respectively. The degenerate unperturbed 
energy is taken to be £^0 = Ij and the interaction V is 
given with the coupling strength x. The relevant matrix 
elements are given by 

and 

Q^Q- 5.)- (49) 




FIG. 3: Dependence of {F„-,{E)} on E with x = 0.2. The 
graphs of y = Fi{E) and y — F2{E) are shown in solid and 
broken lines, respectively. The graph of y = E is also shown. 

We present in Table H] the results of the iterative cal- 
culations, based on Eq. (jiOl) for the KK and Eq. (j42]) for 
the EKK methods. The coupling strength x in the model 
Hamiltonian is taken to be a; = 0.05. Table U shows the 
results for the lowest and second-lowest eigenvalues of H. 
We observe that the convergence in the EKK method is 
much faster than that in the KK method. The EKK 
iteration reaches convergence to fourteen decimal places 
after three iterations. It is also impressive that the KK 
iteration shows linear and steady convergence. As has al- 
ready been discussed in Sec. II, the iteration in the EKK 
method is essentially equivalent to that in the Newton- 
Raphson method which derives quadratic convergence. 
The calculation in the EKK method verifies this the- 
oretical prediction that the number of correct digits is 
doubled at each iteration. 

We depict, in Figs. [2] and |3l the dependence of the 
functions {Gm.iE);m = 1, 2} and {FmiE);m = 1, 2}, re- 
spectively, on the energy variable E with the coupling 
strength x — 0.2. One observes in Fig. [2] that there are 
two poles in the graphs of {Gm{E)} associated with the 
Q-box. On the other hand, the poles disappear in the 
graphs of {Fm{E)} associated with the Z-hox. One sees 
four intersections oiy = Gm{E) and y — E, and six inter- 
sections of y = F,n{E) and y — E ioi m — 1,2. Among 
the six intersections two of these correspond to the pole 
energies, the eigenvalues of QHQ. The energy deriva- 
tives {F'^{E)} at the intersections should be zero for the 
true eigenvalues of H , which are shown in Eq. (j2ip . On 
the contrary the energy derivatives {F/„(_B)} at the in- 
tersections should be two for the solutions of the pole 
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TABLE II: Convergence of the solutions obtained with 
the graphical method for the model Hamiltonian with the 
strength x = 0.05. The E, and F'{Ei) for 1 < i < 4 are the 
true eigenvalues of H and the energy derivatives of {Fm{E)} 
at E = Ei, respectively. The Eq^ and Eg^are the solutions 
for the pole energies. The F^isq-^^ and Fl-^lsq^) are the en- 
ergy derivatives of {Fm{E)} at E = Sq^ and e^^i respectively. 
The notation c indicates convergence to more than six dec- 
imal places. The number of repeats in the modified secant 
method is given by n. The exact values for the solutions are 
also given. 



n 


El 


F'iEi) 


E2 


F'{E2) 


1 


0.890462 


< 10-" 


2.215995 


0.00011 


2 


0.890450 


< 10"" 


2.215681 


< 10"" 


3 


c 


< 10-" 


c 


< lO"'' 


Exact 


0.890450 


0.0 


2.215681 


0.0 


{Ea, Eb) 


(0.8, 0.9) 




(2.2, 2.3) 




n 


E3 


F'iEs) 


Ei 


F'{Ei) 


1 


2.837465 


0.41335 


8.781060 


0.04125 


2 


2.860730 


0.02242 


8.781706 


0.00014 


3 


2.862158 


0.00004 


8.781708 


< 10"" 


4 


2.862161 


< 10-" 


c 


< 10-" 


Exact 


2.862161 


0.0 


8.781708 


0.0 


(Ea, Eb) 


(2.8, 2.9) 




(8.78, 8.80) 




n 




F'iSq,) 




F'iSq,) 


1 


2.750350 


1.98829 


8.751088 


1.95724 


2 


2.749587 


1.99995 


8.750417 


1.99999 


3 


2.749583 


2.00000 


c 


2.00000 


Exact 


2.749583 


2.0 


8.750417 


2.0 


{Ea, Eb) 


(2.70, 2.75) 




(8.70, 8.77) 





energies. The number of the intersections in the graphs 
of {Gm{E)} and {Fm{E)} is verified as already pre- 
dicted in Sec. II. The theoretical predictions on the en- 
ergy derivatives {F^(£')} of these functions are also ver- 
ified as shown in Fig. 3. 

We can obtain, from these graphs in Figs. 2 and 3, 
much information on the convergence conditions in some 
iteration methods. One observes in Fig.[2]that the energy 
derivatives of the first and third low-lying intersections 
are less than one and those of the second and fourth 
are larger than one. From these observations we may 
say that the iteration in Eq. (j40| cannot reproduce the 
second and fourth eigenvalues of H. On the other hand, 
one sees from Fig. [3] that there are four intersections with 
zero energy derivatives. This means that the iteration in 
Eq. (j42|) in the EKK method always converges to any of 
the four true eigenvalues of H. 

We made another numerical calculation by using the 
Z-box and the associated function Fm{E) in the modified 
secant method. Tables |lTl III, and IV show the results 
for the coupling strength x = 0.05,0.1, and 0.2, respec- 
tively. The starting energies Ea and Eb are taken as 
approximate solutions for the intersections of the graphs 
of {y = Fm{E)} and y = E. The results show that aU 
the solutions to Eq. ([TSl) are reproduced. They include 
the four true eigenvalues • • • , £^4} of and the two 
pole energies {£qnSq2\- The values of the energy deriva- 



TABLE III: Convergence of the solutions obtained with 
the graphical method for the model Hamiltonian with the 
strength x = 0.1. As for the notations see Table II. 



n 


El 


F'{E,) 


E2 


F'{E2) 


1 


0.648389 


0.00002 


2.552215 


0.04647 


2 


0.648250 


< 10"" 


2.553840 


0.00016 


3 


c 


< 10"" 


2.553845 


< 10-" 


Exact 


0.648250 


0.0 


2.553845 


0.0 


{Ea, Eb) 


(0.6, 0.7) 




(2.55, 2.65) 




n 


E3 


F'{E3) 


Ei 


F'{E,) 


1 


3.649852 


0.00005 


8.647500 


0.00392 


2 


3.650111 


< 10"" 


8.647794 


< 10"" 


3 


c 


< 10"" 




< 10"" 


Exact 


3.650111 


0.0 


8.647794 


0.0 


{Ea, Eb) 


(3.6, 3.7) 




(8.63, 8.65) 




n 




F'{Sq,) 




F'{eq,) 


1 


2.498646 


1.99121 


8.501897 


1.99693 


2 


2.498334 


1.99999 


8.501666 


2.00000 


3 


c 


2.00000 


c 


2.00000 


Exact 


2.498334 


2.0 


8.501666 


2.0 


{Ea, Eb) 


(2.48, 2.50) 




(8.50, 8.52) 





TABLE IV: Convergence of the solutions obtained with the 
strength x = 0.2. As for the notations see Table II. 



n 


El 


F'{Ei) 


E2 


F'{E2) 


1 


-0.149272 


0.00008 


2.577187 


0.00535 


2 


-0.149586 


< 10-" 


2.579424 


< 10-^ 


3 


c 


< 10"" 


2.579425 


< 10-" 


Exact 


-0.149586 


0.0 


2.579425 


0.0 


{Ea, Eb) 


(-0.2, -0.1) 




(2.5, 2.6) 




n 


E3 


F'{Es) 


E4 


F'{E,) 


1 


5.645253 


0.00003 


8.923645 


0.00244 


2 


5.645051 


< 10-" 


8.925109 


< 10-^ 


3 


c 


< 10-" 


8.925110 


< 10-" 


Exact 


5.645051 


0.0 


8.925110 


0.0 


{Ea, Eb) 


(5.6, 5.7) 




(8.9, 9.0) 




n 


Eqi 


F'{s,i) 


^12 


F'{eq,) 


1 


1.993869 


1.99886 


8.007196 


1.99921 


2 


1.993341 


2.00000 


8.006659 


2.00000 


3 


c 


2.00000 


c 


2.00000 


Exact 


1.993341 


2.0 


8.006659 


2.0 


{Ea, Eb) 


(1.9, 2.0) 




(8.0, 8.1) 





lives {Fl^{E)} at the solutions were also calculated. As 
already proved in Sec. II theoretically, the energy deriva- 
tive at the true eigenvalues of H should be zero and, on 
the other hand, two at the pole energies. These theo- 
retical predictions are also confirmed numerically. From 
this difference of the derivatives of F„i {E) we can easily 
classify the solutions into two parts, the true eigenvalues 
of H and the pole energies. 

As a whole the convergence rates are reasonable in the 
cases with the coupling strength x — 0.05,0.1, and 0.2. 
Three steps of the calculations in the modified secant 
method are enough for yielding the results with accu- 
racy to six decimal places. However, the convergence 
depends strongly on the choice of the starting energies 
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{Ea,Eb). When the spacing of two solutions among {Ek} 
and {Sg} is very narrow, we need to draw graphs accu- 
rately enough for finding approximate values of {Ea,Eb} 
to bracket each of the solutions. 



IV. CONCLUDING REMARKS 

We have introduced a new vertex function called the 
Z-box, which is an operaor defined in a model space and 
a function of an energy E. The Z-box is constructed in 
terms of the Q-box and its energy derivative originated by 
Kuo et al. We have proved that a new expression of the 
effective interaction can be derived by replacing the Q- 
box by the Z-hox in the KK method. With the Z-hox we 
have also introduced a set of scalar functions {Fm{E)}. 
It has been shown that the Z-hox and the associated 
functions {Fm{E)} have the following properties: 

(i) The true eigenvalues {Ek} of the original Hamil- 
tonian H can be given by the roots of a set of 
equations {Fjn{E) — E}. These equations have the 
roots, as additional solutions, at the pole energies 
{Sq} of the Q-box. 

(ii) The Z-hox and the functions {Fm{E)} have no sin- 
gularities even at the poles of the Q-box, and they 
are well-behaved functions of E. 

(iii) The derivatives of {F„i{E)} at the solutions for E 
take two values, i.e., zero at the true eigenvalues 
{Ek} of H, and two at the pole energies {Sq}. 

On the basis of the properties (i), (ii), and (iii), we have 
proposed a new iteration scheme written as = 
FmiE^^^). This iteration always converges and repro- 
duces the true eigenvalues {Ek} of H. Since the energy 
derivatives of {Fra{E)} at the solutions for {Ek} are al- 
ways zero, this method can be understood to be equiva- 
lent to the Newton-Raphson method used to solve non- 
linear equations. The Newton-Raphson iteration derives 
quadratic convergence. Therefore, we can expect that 
this new iteration leads to fast convergence. We have car- 
ried out a test calculation and confirmed the quadratic 
convergence. 



As another method for solving a set of equations 
{Fjn{E) = E}, we have proposed a new non-iterative 
method which we call the graphical method. The solu- 
tions of these equations can be obtained as the energies at 
the intersections of the graphs {y — F„i{E)} and y = E. 
Using the property (iii) we can classify the intersections 
into two parts, i.e., those for the solutions of the true 
eigenvalues {Ek} and of the pole energies {Sq}. These 
graphs make us possible to estimate roughly the posi- 
tions of the roots. With approximate solutions obtained 
from the graphs we proceed to make more accurate calcu- 
lations of the solutions, where we have employed a mod- 
ified secant method which is a combined method of the 
secant and the Newton-Raphson methods. This method 
has been shown to provide a suitable scheme for obtain- 
ing accurate results if we start with approximate energies 
close to the solutions. 

We have made test calculations to assess the graph- 
ical method. We have confirmed numerically that the 
modified secant method reproduces successfully all the 
solutions, including the true eigenvalues {Ek} oi H and 
the pole energies {Sq}. The theoretical predictions in 
the property (iii) on the energy derivatives of {Fm{E)} 
are also verified numerically. We wish to note that the 
graphical method, implemented by the modified secant 
method, yields always convergent results, where we do 
not need any information on the eigenstates and/or the 
eigenvalues of H such as P-space overlaps and/or energy 
spacings. 

We may conclude that the present approach, the 
graphical method with the vertex function Z{E), would 
be promising in resolving some of the difficulties in the 
derivation of the effective interaction. 
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